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Describing and understanding the essence of quantum entanglement and its connection to dynam- 
ical chaos is of great scientific interest. In this work, using information geometric (IG) techniques, 
we investigate the effects of micro-correlations on the evolution of maximal probability paths on 
statistical manifolds induced by systems whose microscopic degrees of freedom are Gaussian dis- 
tributed. We use the statistical manifolds associated with correlated and non-correlated Gaussians 
to model the scattering induced quantum entanglement of two spinless, structureless, non-relativistic 
particles, the latter represented by minimum uncertainty Gaussian wave-packets. Knowing that the 
degree of entanglement is quantified by the purity V of the system, we express the purity for s- 
wave scattering in terms of the micro-correlation coefficient r - a quantity that parameterizes the 
correlated microscopic degrees of freedom of the system; thus establishing a connection between 
entanglement and micro-correlations. Moreover, the correlation coefficient r is readily expressed in 
terms of physical quantities involved in the scattering, the precise form of which is obtained via our 
IG approach. It is found that the entanglement duration can be controlled by the initial momentum 
p , momentum spread a and r. Furthermore, we obtain exact expressions for the IG analogue 
of standard indicators of chaos such as the sectional curvatures, Jacobi field intensities and the 
Lyapunov exponents. We then present an analytical estimate of the information geometric entropy 
(IGE); a suitable measure that quantifies the complexity of geodesic paths on curved manifolds. 
Finally, we present concluding remarks addressing the usefulness of an IG characterization of both 
entanglement and complexity in quantum physics. 

PACS numbers: Probability Theory (02.50.Cw), Riemannian Geometry (02.40.Ky), Complexity (89.70.Eg), 
Entropy (89.70.Cf), Quantum Entanglement (03.65. Ud). 



I. INTRODUCTION 

One of the most debated features of composite quantum mechanical systems is their ability to become entangled 
[HQ- By quantum entanglement we mean quantum correlations among the distinct subsystems of the entire composite 
quantum system. For such correlated quantum systems, it is not possible to specify the quantum state of any subsystem 
independently of the remaining subsystems. 

The generation of quantum entanglement among spatially separated particles requires non-local interactions through 
which quantum correlations are dynamically created 

Quantum entanglement is an indispensable resource for quantum information processes [6|. Continuous Variable 
Quantum Systems (CVQS) are also an interesting topic in quantum information theory By CVQS we refer 
to quantum mechanical systems on which one can - in principle - perform measurements of certain observables 
whose eigenvalue spectrum is continuous. Examples of CVQS are the quantized motion of massive particles with 
the corresponding position and momentum observables and the quantized mode of the electromagnetic field with its 
quadrature observables among others. Most examples of possible applications of entanglement in continuous variable 
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quantum information is based on EPR states [l| or in the optical case, approximations of EPR states using squeezed 
states (il-Hoj. Continuous variable entanglement has been investigated in the context of photon-atom scattering [Tl| . 
photoionization processes [l2l - [l5j . trapped atoms [HI], and classically chaotic systems (l7L Il8j. Correlated CVQS can 
be used as an invaluable non-classical resource for quantum computation and quantum communication Q. 

The most realistic approach to the generation of entangled continuous variable systems is via dynamical interaction, 
of which local scattering events (collisions) are a natural, ubiquitous type [19]. Scattering can result in a decomposition 
of the wave function into transmission and reflection modes. Due to the mutual interactions present in scattering 
processes (such as interference between incoming and reflected parts of the wave function of the composite system), 
quantum particles can become entangled. Moreover, scattering may result in a distortion (due to rapid fluctuation 
of scattering amplitude with relative momentum for instance [3j, |2(J |) of the shape of the wave function. In cases 
with constant amplitudes the wave function of the system may be rendered inseparable as a consequence of reflection 
induced distortion. Interference between incoming and reflected parts of the wave function of the system or the 
distortion effect can result in a non-separable post collision two-particle state. Entanglement generation in non- 
relativistic scattering of distinguishable particles has been investi gate d by a number of researchers 0, [UH^ ■ Most 
treatments consider interactions among similar particle types Il9l - [2ll . l2~il. l25l|. It is however, unclear as to how 
the interaction (scattering) potentials and incident particle energies control the strength of entanglement [3]. As 
will be seen in what follows, the information geometric approach employed in the present work lends some degree of 
clarification on this issue. Furthermore, describing and understanding the complexity of quantum processes is still an 
open problem and our present knowledge on the relations among complexity, chaoticity and quantum entanglement 
are not at all satisfactory 27]. As we will see, our work sheds some lights on this issue as well. 

In this article, we explore the potential utility of the Information Geometric Approach to Chaos (IGAC) [28 
|33| for analyzing quantum mechanical systems. The IGAC is a theoretical framework developed to study chaos in 
informational geodesic flows on statistical manifolds associated with probabilistic descriptions of physical systems. 

We seek to provide a quantitative estimate of the degree of entanglement of CVQS in terms of information geometric 
quantities such as solutions to geodesic equations (expected values of momentum and momentum spread in this case) 
and micro-correlation coefficient r. The quantity r parameterizes the correlated microscopic degrees of freedom of the 
system. An important question that arises is whether or not r can be understood in terms of physically measurable 
quantities. 

As described above, when particles with no initial correlations collide, they may emerge from the interaction 
entangled [34]. Hence, we consider two CVQS with Gaussian continuous degrees of freedom that are prepared 
independently, interact via a scattering process mediated by an interaction (scattering) potential with finite range 
and separate again. We investigate the entanglement of the two-particle wave function of the system generated by 
such a scattering event. In this context we ask the question: how much entanglement between the two particles is 
generated and on what does it depend? Surprisingly, there are only a few studies of entanglement production from 
the scattering of two particles (35L l36| . The nature of quantum entanglement arising from s-wave scattering has yet 
to be fully explored [37j • We choose to use Gaussian states (H, [38| since many important properties of these states 
can often be obtained in an analytic fashion. Moreover, it is known that a good way to describe naturally occurring 
quantum states is as spatially localized Gaussian wave-packets, or as density matrices built from them (23 . I39TI421 ] . 

For a system of two spinless, structureless, non-relativistic particles with no internal degrees of freedom, a complete 
set of commuting observables is furnished by the momentum operators of each particle [34j . The continuous variables 
in our case are therefore taken to be the momentum of each particle. We investigate how the initial conditions and 
the magnitude of r of the system affect entanglement; specifically, the duration of entanglement. By duration of 
entanglement we mean the temporal behavior of the magnitude of entanglement. The IGAC is also used to estimate 
the extent to which the complexity of the geodesic information flow (continuous spectrum of expected values of some 
relevant observable) of the quantum system is affected by r. 

The layout of this article is as follows. In Section UH we reexamine the s-wave scattering induced quantum en- 
tanglement of two spinless, structureless, non-relativistic particles, which are represented by two-particle Gaussian 
wave-packets (24|. We exploit the fact that the three-dimensional scattering problem can be effectively reduced to 
a one-dimensional problem as far as the three-dimensional representations of wave-packets are isotropic. In Section 
IIII1 we outline the main ideas behind the IGAC and present the information geometry of correlated and uncorrelated 
Gaussian statistical manifolds, which is to be employed in our investigation of quantum entanglement and complexity 
in the following sections. In Section ITVl we use information geometric techniques in conjunction with standard partial 
wave quantum scattering theory to provide an information geometric characterization of quantum entanglement. In 
Section |V] we obtain exact expressions for the information geometric analogue of standard indicators of chaos such 
as sectional curvatures, Jacobi field intensities and Lyapunov exponents. Finally, we present an analytical estimate 
of the information geometric entropy (IGE) and this allows us to connect quantum entanglement to the complexity 
of informational geodesic flows in a quantitative manner. Our concluding remarks are presented in Section [VII 
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II. ENTANGLEMENT AND GAUSSIAN WAVE-PACKET SCATTERING PROCESSES 



In this Section, we reexamine the s-wave scattering induced quantum entanglement of two spinless, structureless, 
non-relativistic particles, represented by two-particle Gaussian wave-packets as presented in [24]] ■ For ease of analysis, 
we exploit the fact that the three-dimensional scattering problem can be effectively reduced to a one-dimensional 
problem as the three-dimensional representations of wave-packets are isotropic. 



A. The Pre-collisional Scenario and the Effective Dimensional Reduction 



For the purpose of modeling a head-on collision we consider two identical (but distinguishable), spinless particles 
in momentum space, each represented by minimum uncertainty Gaussian wave-packets. Before collision, particles 1 
and 2 are initially located far from each other - a linear distance R - each having the initial average momentum 
(pi) o = p and (p2) = — Po, respectively, with equal momentum dispersion a a (see Figure [T]). The normalized, 
separable (i.e., non-entangled) two-particle Gaussian wave function representing the situation before collision is then 
given by [24j 



V>(ki,k 2 ) = Vi (ki)®V 2 (k 2 ), 
with respective single particle wave functions 

^ 1/2 (k 1/2 ) = a (k 1/2 , <k 1/2 ) o ; a ko ) e K^~M )^ 

where 



a(k 1/2 ,(k 1/2 ) o ;(T fco ) = 



/ 1 



3/4 



exp 



(k 1/2 - (k 1/2 ) Q )' 

Mo 



(1) 

(2) 
(3) 



with k 1/2 = ^i, (k 1/2 ) o 



±2^- = ±k Q , ako = ail d qi/ 2 = T^R-o- Observe that wave functions in @ 

(4) 



satisfy the normalization conditions 

J (kO^kOd 3 ^ = J 'a 2 (k^ki^fffaO^ki = 1 



and 



J i>* 2 (k 2 ) V 2 (k 2 ) d 3 k 2 = J a 2 (k 2 , (k 2 ) Q ; a ko ) d 3 k 2 = 1. 



(5) 



The type of state described by (p} is ubiquitous when describing quantum systems of continuous variables. 

One should note that the three-dimensional Gaussian wave-packet (fTJ) is isotropic. That is to say, in polar coordinates 
the representation of the separable two-particle state exhibits a functional dependence on the radial variable only. 
For this reason the three-dimensional vectorial representation (fT]) may be effectively reduced to a one-dimensional 
representation. This may be demonstrated as follows. First, using (fTJ), ([2]) and (J3j) , we express the two-particle 
wave-packet in Cartesian coordinates as 



ip , kf, kf ) , (k 2 , k 2 )) — ipi (k 1 ,k\ 1 k 1 ) (g) ij) 2 (k 2 , fc 2 , k 2 ) , 



where 



^1/2 (^l/2> ^1/2' ^1/2) 



2™L 



3/4 



n exp 



( k i/2 ( fc i/ 2 ) o 



= *( fc l/2 ( fc l/2) o )9l/2 



(6) 



(7) 



with fc^y 2 = kx/2 • e 



u 

J' \ 1/2 



= (k!/ 2 ) o • e^, q^y 2 — qi/ 2 • e 3 , the components in a Cartesian basis {e^}. Then the 
probability density for this wave-packet becomes 

iv (k 1; k 2 )i 2 = r m, kf, kD , (fcf , ki m $ m,k\, ^) , (^, ki m 



n exp 



k{ - (ki 



2< 



exp 



u _ I u 
fi, 2 i h, 2 



2< 



(8) 
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Figure 1: An illustration of the two-particle system before and after a head-on collision. Before collision the two particles are 
initially the distance R away from each other and move toward each other with the momenta p and — p , respectively; both 
particles have the identical momentum spread a (h/a in configuration space). After collision a large spherical shell represents 
the scattered part of the single particle density (i.e. either particle 1 or particle 2) under the s-wave approximation. The arrows 
indicate that after collision the two particles move away from each other with the momenta, — p Q and p Q , respectively 



Upon integrating \ip (k i; k2)| 2 over d 3 kid 3 k2, one obtains 

+oo +oo -t-oo -t-oo +oo +oo 



////// II dM^hMkl,k2)| 2 = /l/2 
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3 =X,y,z 
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where 



2™L 



3/2 



-oo +oo + 00 



dkj exp 
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with 1 = 1,2. 



(10) 



Converting the integrals // (Z = 1,2) into polar-coordinate representation yields 



2^L 



2^L 



3/2 



-oo +oo +oo 



3/2 
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+00 ^ 

47r / dkikf exp ( - 1 
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(11) 



where kf = = (kf ) 2 + (kf) + (kf ) 2 (I = 1, 2). The remaining integral in (jTTJ) can be replaced with 

+00 



-I-00 +00 

/ dkikf exp {Si) = ^ / * cxp (-4:) 





+00 



— — \ dk\ exp . 
2 / P \ 2<t 2 



ko 



+00 



/ dfc ; exp 



(fcj - (fcQo) 2 
2aL 



(12) 



where we have changed the domain of fc; from [0, +00) to (—00, +00) in order to obtain the last equality. By 
substituting (fl~2|) into (fTTjl . followed by substituting (fTTj) into the right-hand side of (j9|) and finally inserting ([8]) into 
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the left-hand side of we establish 
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The integration of |^/> (ki,k2)| 2 over <i 3 ki<i 3 k 2 now reads 



|^(ki,k 2 )| 2 d 3 kid 3 k 2 



2™L 



exp 



(ki - (kQJ 2 + (k 2 - (k 2 )J 2 
2aL 



d 3 kid 3 k 2 



+ 00 + 00 



2™L 



exp 



— OO —CO 



(fci - (fcQJ 2 + (fc 2 - (fc 2 ) o ) 2 
2a 2 feo 



dfcidfc 2 . (14) 



Thus, we may effectively reduce the three dimensional two-particle wave function if) (ki, k 2 ) expressed via (JT]), ([2|) and 
([3]) to the two-particle one-dimensional wave function ip (fci, fc 2 ) given by 



V» (fci, fe) = -01 (fci) ® ^2 (fo) , 
with respective single-particle wave functions, 

0i/ 2 (fci/ 2 ) = a (fci/2, (fci/ 2 > ; °ko) e'^-M.)^ 

where 



a(fci/ 2 ,(fci/2> o ;^o) = (^~2~) 



1/4 



exp 



(fcl/2 ~ (fcl/ 2 ) Q ) 

44o 



(15) 



(16) 



(17) 



and ki/2 = € (— oo,+oo), (fei/a) = V = =^^f = i&oj 0>o = X' ^1/2 = ^l-^o- The wave functions (flU)) 
satisfy the normalization conditions 



+ CO 



-Hoc 



V'l/a ( fc l/2) 01/2 (fcl/2) *l/2 = / a 2 (fci/2, (fcl/ 2 ) ; CTfco) dfci/2 = 1- 



(18) 



B. The Post-collisional Scenario 



After collision, the wave function for the two-particle system in the long time limit takes the form |24| : 

V (ki, k 2 , t) = (Ny 1/2 k (ko ^ 2 (k 2 ) e-^K^+^Dvprn) + e?/ , scat (ki) k2) t y 



(19) 



where N and £ are normalization constants such that tp and V'scat are both normalized to unity, m denotes the mass 
of each particle. Here t/'i/a (^1/2) i s given by ([5]) and ((31) and thus we have 



il> x (ki) ?/> 2 (k 3 ) e -«( fc ?+ fc 2) t /( 2m ) 



/ ! X3/2 



(kr - k G ) 2 + (k 2 + k G ) 2 



Ho 



(20) 



xe -i{k 1 -k a )R a /2+i(k 2 +k a )R /2-ih{k 2 1 +kl)t/{2m) 

We treat |e| -C 1 as a small number. Following [24j|, one can write 

e^ Bcat (k 1 ,k 2 ,t)=eV c . m .(K,t) >7 (k,t), (21) 
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where 



V>c.m. (K,t) = 

and the scattering part erj (k, t) is given by 

1 



/ 1 



3/4 



exp 



K 2 

8al 



-ihK 2 t/{2M) 



erj (k, t) 



n ik' r 



V> re ] (k',0) / (fc') e -*-r- l r i fc' 2 t/(2 M ) d 3 k / d 3 r 



with 



VVel (k',0) = 



3/4 



exp 



(k' - kj 



,-%{k-ko)Ro 



(22) 



(23) 



(24) 



Here we have adopted the center of mass and relative coordinates such that the conjugate momenta K = ki + k2 
and k = i (ki — k 2 ) are used along with the total mass, M = 2m and the reduced mass /i = m/2. The quantity 

/ (k) = - — ^il^ 1 ^ s ^ ne s_wave scattering amplitude due to the s-wave scattering phase shift 9 (k). Inserting (|24p into 
and performing the integral, we obtain 



£77 (k, t) 



3/4 



exp 



2-L 



g (k) e 



-i(k-k a )R -ihk 2 t/(2^) 



where 



4» (k -ial R ) fc 2 /(fc) 

?(*) = 2 • 



(25) 



(26) 



and the approximation has been made under the assumption of low energy s-wave scattering. Then by (|21[) . (1221) and 
(j2"5) we find 



eV'scat (ki,k 2 ,i) 



2™L 



3/2 



exp 



K 2 + 4 k - 



(k-fc k)' 



X £ (fc) e -i( k - k o)Ra-ihK 2 t/(2M)-thk 2 t/(2ti) 



Due to the fact that 



exp 



exp 



(k x - k Q ) 2 + (k 2 + k D r 



4< 



-i(k 1 -k )R /2+i(k 2 + k )R /2-ih(kl+kl)t/(2m) 



K 2 + 4 k - fc n k 



~\ 2 



J -i(k-k )R -ihK 2 t/(2M)-ihk 2 t/(2 l i) 



we may combine (j2"0")) and (f2"Tf to rewrite (fH?)) as 

^(kx.ka.t) = (Ny 1/2 



2™L 



3/2 



exp 



K 2 +4 



(k-fc G k) 2 



8< 



X [1 + gf k ^ e -i{k-k )R -ihK 2 t/{2M)-ihk 2 t/(2fj,)_ 



(27) 



(28) 



(29) 



Our dimensional analysis carried out in the separable case (see the previous Subsection) applies equally well in the 
entangled case. Hence, we may reduce the three-dimensional wave-packet rp (ki, k 2 , t) expressed via (fT9"|) to the 
one-dimensional one, 



(30) 
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Pl/2 



where the single-particle wave function ip 1 / 2 (^1/2) is specified via (JTHJ) and (|T7)) . together with £4/2 — — ^~ 
(-00, +00), (h/ 2 ) = ^ Pl/ f^° = ±x = ±fco, cr fco = g x/2 = =f|-Ro- In analogy to (J^SJl, i/)(ki,k 2 ,t) can be 



rewritten as 



i>(h,k 2 ,t) = (N) 



-1/2 



2™L 



1/2 



exp 



K 2 + 4 (fc - fc Q )' 
8< 



X [f + Q(k)] e - i ( fc - fc °) fl =- iriK2 */(2M)-ttfc 2 t/(2A') j 



(31) 



where we adopt the one-dimensional center of mass and relative coordinates, whose conjugate momenta are defined 
as K = k\ + /c2 G (—00, +00) and A: = | (fci — fc^) G (— 00, +00), and p(fc) is given by f|26[) . Separating variables in 
T]), we may write 



iv(fci,fc 2) i)r 



il>(ki,k 2 ,t) ip* (ki,k 2 ,t) 
K 2 



N- 1 

2— exp 



27TCT 



/CO 



exp 



i + 2»( e (fc)) + | e (fc)r 



(32) 



where fc = fc — fc and -ft denotes the real part of p (fc). However, we find that the complex- valued scattering amplitude 
/ (k) can be approximated as real since / (k) = + O (# 2 ) for S (fc) < 1. In view of this fact and (|2"o| we employ 
the following approximations: 



$t(g(k)) « 4i? fc 2 /(fc), 



(33) 



16(fc 2 + a|X) fc 4 (/W) 2 



(34) 



feo 



Upon Taylor expanding / (k) and (/ (A:)) 2 around k — fc Q , we re-express 3ft (p (k)) in (f33|) and |p (A;)| 2 in (|34|) as 

00 

H (e (fc)) « 4i? Q fc 2 A n (k - kof (35) 

n=0 

and 

| ff (fc)i 2 « 16 ( fc ° + fc4 f; a. (* - fc ) n , (36) 



'feo 



ri=0 



respectively. Here the quantities A n and -B„ are appropriate coefficients determined from the expansions / (fc) 
S^=o A n (k - k Q ) n and (/(fc)) = J2n=o B n {k - feo)™, respectively. By inserting (J3SJ) and (f3SJ) into 
integrate |^| 2 as follows: 

+00 +00 +00 +00 



\ip\ 2 dkxdk 2 = 2na 2 ko N / / |^| 2 dfTdjfc 



— OO —OO 



— OO —OO 



dK exp 



"4a 2 D 



a/cexp — 



ko 



00 

1 + 8i? (fc 2 + 2fc Q fc + fc 2 ) A n ~k r ' 



n=0 



16 (fc 2 + oi R 2 ) 



OO 

(fc 4 + 4fc fc 3 + 6fc 2 fc 2 + 4fc3fc + fc 4 ) B n k r ' 



71=0 



(37) 



where k = fc — fc Q . From [43j we find 

+00 

c?fc exp 



fc" = (2m - l)\\^a ko 



'ko 



(38) 



where m = 0,1,2,.... It should be noted however, that we have approximated / (fc) as a real-valued function 
/ (fc) f=a assuming 9 (fc) <C 1. For low energy s-wave scattering, which is the case presently under consideration, 
we have k <C 1 and 6 (k) = —ka s + O (fc 2 ), where the parameter a s of dimension length is defined as the s-wave 
scattering length (44|. This leads to / (fc) — a s and (fc) 0, where the superscript ^ denotes any p-th. order 
derivative (p = 1, 2, . . .). Hence, we have Ao = — a s , = A2 = ■ ■ ■ = and Bg = a 2 , B\ = _B 2 = ■ ■ ■ = 0. Making use 
of these coefficients as well as (|38|) , one may compute the integral in (|5T|) to obtain 



+00 +00 



+00 +00 



|V>| 2 *idfc 2 = 2 7 r<7 2 . JV 



\ip\ 2 dKdk 



1 4 ^ , 2 s R ^ 4 (fc 2 + ^gg) (4fc Q 4 + 12fc 2 q 2 Q + 3a 4 fco ) 2 
1 - 4 (2fc Q + CT fco ) i? D a s + — ^ ^ — 2 '-at 



ko 



.(39) 



Indeed, one would obtain the same result as (|39[) to the leading order by evaluating the following integral: 

+00 +00 



dk\dk2 exp 



1 



2 (l-^Q M ) I 



(fci - fc D ) 2 (fcx - fc Q ) (fc 2 + fc Q ) (fc 2 + K) 2 

J2 2r QM ~2 1 ~2 



'ko 



2iraLJl-r 2 QM = 2<kq\ 



QM 



Or, 



QM) 



(40) 



where 



r QM = ^8 (2fc 2 + cr 2 D ) i? a s < 1. 
Therefore, to a good approximation, we may replace the probability density in (|32[) with 



(41) 



after 

QM = 



^after ^ fc^j) 



exp 



2 (!-^m) 



(fcl-fco) 2 

— ^ — 



^feoVi-^QM 



(42) 



where the integral in (|40| has been normalized. In case tqm = 0, (|42|) reduces to 



r>after pbefore 



^^"(fcl,^) 



2™L 



exp 



(fc! - fc f + (fc 2 + k Q y 

2a 2 Q 



(43) 



which is verified via §U$, (jTHJ) and (fTTj) 

The obtained expressions for the prol 
tigation as will be explained in the next Section 



The obtained expressions for the probability densities -Pqm™ 6 and Pqm 1 motivate our information geometric inves- 



III. THE INFORMATION GEOMETRIC PERSPECTIVE 



In this Section, we outline the main ideas behind the IGAC and present the information geometry of correlated and 
uncorrelated Gaussian statistical manifolds employed in our investigation of scattering induced quantum entanglement. 



A. On the IGAC 

IGAC [45|, |4(| is a theoretical framework developed to study the complexity of informational geodesic flows 
describing physical. It is the information geometric analogue of conventional geometrodynamical approaches to chaos 
[291 . [30L lr7H49j | where the classical configuration space is replaced by a curved statistical manifold with the additional 
possibility of considering chaotic dynamics arising from non-conformally flat metrics. Additionally, it is an information 
geometric extension of the Jacobi geometrodynamics (the geometrization of a Hamiltonian system by transforming it 
to a geodesic flow [50]). 
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More specifically, IGAC is the application of entropic dynamics (ED) [5l| to complex systems of arbitrary nature. 
ED is a theoretical framework that arises from the combination of inductive inference (Maximum Relative Entropy 
methods, [EH, [EH]) and Information Geometry (IG), that is, Riemannian geometry applied to probability theory 
[54| . IGAC extends the applicability of ED to temporally-complex (chaotic) dynamical systems on curved statistical 
manifolds and relevant measures of chaoticity of such an IGAC have been identified [45j . 

The essential ideas underlying the IGAC and the construction of statistical manifolds are presented in what follows. 
Let the probability distribution function (PDF) P (X\Q) represent the maximally probable description of the system 
being considered. The quantity X is a random variable that represents a microstate of the system, while represents 
a macrostate. The sets {X} and {0} form the microspace X and the parameter space T>q, respectively. The set 
of probability distributions forms the statistical manifold A4 . A geodesic curve on a curved statistical manifold M. 
represents the maximum probability path a complex dynamical system explores in its evolution between initial and final 
macrostates Oi and Op, respectively. Each point of the geodesic on an n-dimensional statistical manifold M represents 
a macrostate parametrized by the macroscopic dynamical variables {0}. Furthermore, each macrostate is in a one-to- 
one correspondence with the probability distribution P (X\&) representing the maximally probable description of the 
system being considered. The main goal of an ED model is that of inferring "macroscopic predictions" in the absence 
of detailed knowledge of the microscopic nature of the arbitrary complex systems being considered. More explicitly, 
by "macroscopic prediction" we mean knowledge of the statistical parameters (expectation values) of the probability 
distribution function that best reflects what is known about the system. This is an important conceptual point. The 
probability distribution reflects the system in general, not the microstates. Once the microstates have been defined, 
we then select the relevant information about the system. In other words, we have to select the macrospace of the 
system. We emphasize that knowledge of both initial and final macrostates is not necessary to carry out macroscopic 
predictions. For instance, one may only have knowledge of the initial state and assume that the system evolves to 
other states, without actually knowing what the final state is. In such a case, it can be shown that the system moves 
continuously and irreversibly along the entropy gradient [55j. We note that in its present form the IGAC can only be 
applied to CVQS. This restricted applicability is due to the fact that the IGAC is used to understand the evolution 
of continuous trajectories on M.. In the context of quantum mechanical systems, the set {0} would correspond to 
the continuous eigenvalue spectrum of some observable (expected values). The IGAC must be reformulated in order 
to be applicable to general quantum systems. Such a reformulation is currently in progress. 

For a brief overview of some of the latest applications of the IGAC to both classical and quantum scenarios, we 
refer to (5rj |. 



B. Gaussian Statistical Models and Micro-correlations 



Here, we introduce the notion of Gaussian statistical models (manifolds) in either absence or presence of correlations 
between the microscopic degrees of freedom of the system (i.e. micro-correlations). 



1. Statistical Models in Absence of Micro-correlations 

Consider a statistical model whose microstates span a n-dimensional space labeled by the variables {X} = 
{xi, X2, ■ ■ ■ , x n } with Xj G R, Vj = l,...,n. We assume the only testable information pertaining to the quanti- 



ties Xj consists of the expectation values 



and the variance Axj. The set of these expected values define the 



2n-dimensional space of macrostates of the system. Our 2n-dimensional statistical model represents a macroscopic 
(i.e. probabilistic) description of a microscopic, n-dimensional physical system evolving over a n-dimensional (micro) 
space. We assume that all information relevant to the dynamical evolution of the system is contained in the probability 
distributions. For this reason, no other information is required. Each macrostate may be thought as a point of a 
2n-dimensional statistical manifold with coordinates given by the numerical values of the expectations and 
The available relevant information can be written in the form of the following 2n information constraint equations: 



+00 



(1) V 2) *j 



Axi 



+00 



dx 



{x j -{x j )fP j (x j \^ 2 H j 



(44) 



The probability distributions Pj {pj^'&j-, $j) in (O are constrained by the conditions of normalization, 

+00 

dXjP^X^^^U^l. 



(45) 
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Maximum Relative Entropy methods [52, [5J, [57], [58| allow us to associate a probability distribution P (X\@) to each 
point in the space ol states {6}. The distribution that best reflects the information contained in the prior distribution 
m (X) updated by the information (ixj) , Axj) is obtained by maximizing the relative entropy 



5(9) = - / dXP(X\e)\n 



P(x\e) \ 

m(X) J 



(46) 



As a working hypothesis, the prior m (X) is set to be uniform since we assume the lack of prior available information 
about the system (postulate of equal a priori probabilities). Information theory identifies the Gaussian distribution 
as the maximum entropy distribution if only the expectation value and the variance are known [59]. Indeed, upon 
maximizing (14"6"|) given the constraints (14"4"|) and (|4"S")) . we obtain 



n 

P(X\Q) = I[P j (x j \W* J W4 j ), 



where 



Pi 



exp 



(47) 



2a 2 



(48) 



and in standard notation for Gaussians, ^fij = (xj) = fj,j, ^fij = Axj = o~j. The probability distribution (|4"Tj) 
encodes the available information concerning the system. 

The statistical manifold Ai associated to (14"T1) is formally defined as follows: 



M = { P [X\Q) = HPj {,-, llr a ,) : 9 = . . 
j'=i 

The parameter space 2?^ otal ) (homeomorphic to A4) is defined as 



V, 



(total) def 



2n 



)T, d k = 



k=l 



,tf 2n ) e2?i total) 



P(X\Q) > 



(49) 



(50) 



where I$k is a subset of R and represents the entire range of accessible values for the macrovariable <d 
The line element ds 2 arising from (j47|) is [60] 



ds 2 M = g ab (9) dTdd b 




with a.b = 1, 



,2n. 



(51) 



Note that we have assumed uncoupled constraints among microvariables Xj. In other words, we assumed that infor- 
mation about correlations between the microvariables need not to be tracked. This assumption leads to the simplified 
product rule (|4"T1) . 

A measure of distinguishability among the macrostates of the Gaussian model is achieved by assigning a probability 

distribution P (X\Q) to each 2n-dimensional macrostate 9 = f { S 2 ^ $j) } n irs == {(i x j) > ^- x j)} n pairs' 
process of assigning a probability distribution to each state provides M. with a metric structure. Specifically, the 
Fisher-Rao information metric g a b (9) [HH is a measure of distinguishability among macrostates on the statistical 
manifold Ai, 



g ab (9)- J dXP(X\0)d a lnP{X\e)dblnP{X\e) -4 y dXd a ^P {X\Q)d b ^JP (A|9), 



(52) 



with a, b = l,...,2n and <9 a = ofa- It assigns an information geometry to the space of states. The information 
metric g a t (9) is a symmetric and positive definite Riemannian metric. For the sake of completeness and in view of 
its potential relevance in the study of correlations, we point out that the Fisher-Rao metric satisfies the following 
two properties: 1) invariance under (invertible) transformations of microvariables {X} £ X\ 2) covariance under 
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reparametrization of the statistical macrospace {0} € T>q. The invariance of g a i, (0) under reparametrization of the 
microspace X implies that [54|, 



X C R" 3 x i— > v d = / (as) e y C M n => p (a|<?) i— > p' = 

x=f-Hv) 

The covariance under reparametrization of the parameter space £>e (homeomorphic to A4) implies that [54 



1 



Pe9^/ d =i f / (0) e Pey =^> to (0) — > <4 (*') = 



9cd {#) 



where 



9'ab K) = / dxp' (x\0') d' a lnp' (x|tf') ^ lnp' , 
with 9; = ^ and p' (asji?') = p (x|tf = r 1 (<?'))• 



(53) 



(54) 



(55) 



2. Statistical Models in Presence of Micro-correlations 

Coupled constraints would lead to a "generalized" product rule in (IT71) and to a metric tensor (|5"!2)) with non-trivial 
off-diagonal elements (covariance terms). In presence of correlated degrees of freedom {xj}, the "generalized" product 
rule becomes 



Ptotal (X 



1 ) • ■ • ) x n 



3=1 



correlations 



' Ptotal ( x l, • • • > x n) £ II Pj(Xj) , 

3=1 



(56) 



where 

Ptotal =-Pf» (^nl^lj • • -,Xn-l) Pn-1 (a: n _i|xi, . . .,X n - 2 ) ■ ■ ■ P 2 (tC 2 |xi) P\ (xx) . (57) 

For instance, correlations in the degrees of freedom may be introduced in terms of the following information-constraints, 

x j = fj • • • ) x j-i) > Vj = 2, . . . ,n. (58) 

In such a case, we obtain 

Ptotal ( x i,- ■ ■ ,x n ) = 5 (x n - f n (xx, . . .,X n -i)) S (x„_i - /„_! (xx, ■ ■ ■ ,X n - 2 )) ■ ■ ■ 6 (x 2 - f 2 (xx)) Pi (xx) , (59) 
where the j-th probability distribution Pj (xj) is given by 



(60) 



Correlations between the microscopic degrees of freedom of the system {xj} (micro-correlations) are conventionally 

(m 

def \XiXjl {Xi} (Xj) 



introduced by means of the correlation coefficients r^ lc 16 1H 



(micro) 



r {xi , Xj ) 



with o i = J ( (xi - (xi)) ), 



(61) 



with 



(micro) 



6 (— 1, 1) and i,j — 1, . . . ,n. For the 2n-dimensional Gaussian statistical model in presence of micro- 



correlations, the system is described by the following probability distribution P (X|0): 



P(X\Q) = 



1 



[(27r)™detC(9)] ! 



■ exp 



-- (X - Mf ■ C' 1 (6) • (X - M) 



3=1 



2oJ 



, (62) 



where X = (xx, • ■ • , x n ), M = . . . , /i„) and C (0) is the (2n x 2n)-dimensional (non-singular) covariance matrix. 

In what follows, we will introduce the three-dimensional micro-correlated Gaussian statistical model being investi- 
gated. 
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C. The Two- variable Micro-correlated Gaussian Statistical Model 



Consider micro-correlated Gaussian statistical models with 2n = 4. For n = 2, (|62f leads to the probability 
distribution P (x, y\fJ, x , fJ, y , o x , er y ) (6lT |. 



P (x,y\fJ. x ,H y ,CT x ,(Ty;r) 



ex Pi ~3p=py 



27TCT X CT V V / 1 — 



bivariate normal distribution where > 0, <r y > 0, r G (— 1, 1). X = (x.y), Q = (n x , f/> y ,cr x ,cr y ) and C (O), 



with f, j = 1,2. 



Substituting (|63l) in (|52[). the Fisher-Rao information metric g a b (n x , fi y ,(T x ,<jy;r) becomes 



(63) 



(64) 



/ i 







_„ : : -" 







<rj(r 2 -l) 





1 



3ab {^ x ^y^x,O y \r) 

The infinitesimal line element ds^ /[4D relative to g a b (/U x , [i y) cr X) <T y ; r) is given by 



o- x (T H (r 2 -l) 





fl 2-r 



(65) 



ds Al*?„ = 3ii {<J X ; r) dfxl + 333 [cr y ;r) d^ y + 322 (<J x ;r) da 2 x + 344 (ov; r) dcr 2 + 2313 (er^, a,,; r) d/i^d/^ 
+2324 (c* , <Tj,; r) da x da y , 



(66) 



where 



5ii (ex!? - ) 



<7|(r»-l) 



#13 (0"a;,0"w; r ) 



VxVy (H - 1) 



2-r 2 



a* ( r 2_l)' 
1 



.924 



{(Jx ^ r) = ^ (r* - 1) ' 331 (a x ,a y ;r) = — p - 1} , 333 (a„;r) = p - 1} 



542 (a x ,a y ;r) 



OxOy (r 2 - 1) 



, 344 (cr y ;r) 



2-r 2 



al (r 2 -l)' 



(67) 



It is rather difficult to present an analytical study of the IGAC associated with infinitesimal line element ds 2 Mm in 
(|66p . Such a study will be the subject of forthcoming investigations. In the present work we consider a special class 
of Gaussian models, namely those in which a y — a x — a. Then the probability distribution P (x, y\^ x , H y , cr x , a y \ r) 
in (l63|) can be reduced to a simpler form, 



CX Pi -2(1^) 



where a > 0, X = (x,y), 8 = (fx x ,fi y ,a) and C (0), 



Ci 



^ Zr ^ 1 ^ 



271-aVl - r 2 



ro 1 + n x n 



with i, j = 1, 2. 



The Fisher- Rao matrix g a b (ji x , (i y ,cr; r) associated with P (x, y|/i x , /x , <j; r) reads 



(68) 



(69) 



3ab (n x ,H v ,cr;r) = \ 



"y, 

^ . 



r 2 — 1 r^" 1 

4 



(70) 
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The line element associated with metric g ab (^ x , fi v , o~; r) is given by 



ds 2 MaD = gn(vx;r)dfil+g 33 (a y ;r)diil+2g 13 (a;r)d^ x dv y + [g22{<J--,r)+g4 4 (a;r) + 2g24(<J;r)}d(T 2 (71) 



1 



1 



rdfil 



1,2 2r 



-dp, x d\x y + 4da' 



a 2 y 1 _ r 2 m ' 1 _ r 2 i_ r 2 
Observe that in the absence of micro-correlations, the two-variable probability distribution (|68p reduces to 



1 



2tto- 



■ exp 



2a 2 



exp 



2a 2 



(72) 



while the metric and corresponding line element become 



1 



1 



and 



9ab (n x ,fJ- y ,a) = — 10 
4 



(73) 



(74) 



respectively. In what follows we limit our analysis to the study of non-negative micro-correlations, that is, we will 
consider re [0, 1). 

D. The Information Dynamics on the Statistical Manifold M^orr. 

The information dynamics on the manifold At^Sr represented by (|70p can be derived from a standard principle of 
least action of Jacobi type [5l|. The geodesic equations for the macrovariables of the Gaussian ED model are given 
by nonlinear second order coupled ordinary differential equations, 



H^ 2 



„ dd b dd c 



dr dr 



(75) 



where a,b,c= 1, 2, 3 and we denote $ = ^ = fj, x , d = /i 2 = l^ y , 'd — <J. The connection coefficients T^ c appearing 
in (1751) are defined as 1621 



(76) 



r? c = \g ad {d b g dc + d c g bd - d d g bc ) . 

In our case, through (|T0|) the non- vanishing connection coefficients are given by 

1 _„ r 



T^l " T~l2 ^" T~l2 

1 13 - 1 23 - - 1 32, 1 11 - - 4(T (r2 _ 1} 



r 3 - 

) 1 12 — 



4cr (r 2 - 1) 



T^l3 T^3 

— 1 21, 1 22 — ~ 



1 3 _ 1 

4a(r 2 -l)' ^''a' (7?) 



The geodesic equations in (|75[) describe a reversible dynamics whose solution is the trajectory between an initial 
0i and a final macrostate Of- The trajectory can be equally well traversed in both directions. In the case under 
consideration, substituting (|77[) in (|75[) . the three geodesic equations become 



d 2 ^ (r) 


2 


d^, 1 (r) da (r) 


dr 2 


a(r) 


dr dr 


d 2 n 2 (r) 


2 


dfi 2 (r) da (r) 


dr 2 


a(r) 




d 2 o (r) 




/da(r)\ 2 


dr 2 


a(r) 


\ dr J 4 



d^i ( T ) 



d^2 ( T ) 
dr 



(78) 
(79) 
(80) 



d/i-t (r) d/x 2 (r) 



2cr(r)(r 2 -l) dr 



dr 
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Integration of the above coupled system of nonlinear differential equations is non-trivial. A detailed derivation of 
the geodesic paths is given in Appendix [A"l After integration of (|78p . (|?9")) and (f50j) . the geodesic trajectories for the 
non-correlated Gaussian system become, 



!(t;0) = -Vpf+2o|tanh(A r), 
/i 2 (r;0) = v/p2 + 2a2tanh (A Q r) , 

a(T;0) = cosher) ' 



while for the correlated Gaussian system the geodesies read, 



where the subscript " " denotes the initial state, and 

A = — sinh -1 [ — — 
To V V2cr 



Ml (r;r) = - y/(l-r) {pi + 2ag ) tanh (A Q r) 
/i 2 (r;r) = - r) (p§ + 2^) tanh {A oT ) , 

TZ I i 



2 P ° + <J ° cosh (A oT )' 




whose derivation is presented in Appendix [Bl 



(81) 
(82) 

(83) 



(84) 
(85) 

(86) 



(87) 



IV. APPLICATION OF INFORMATION GEOMETRY TO QUANTUM PHYSICS - PURITY, 
SCATTERING AND QUANTUM ENTANGLEMENT 

In this Section we use information geometric techniques in conjunction with standard partial wave quantum scat- 
tering theory to provide an information geometric characterization of quantum entanglement. 



A. Association of Quantum Systems with Information Geometric Systems 



We now focus on applying IG methods to the quantum entanglement produced by a head-on collision between 
two Gaussian wave-packets in momentum space. We observe from (| 1 5 1) and ([T^)l that the two-particle probability 

2 

density before collision is given by ip e ° re (fcb fe) = a 2 {ki, (ki) D ; <7fc ) a 2 (&2, (^2)0 I Cfc ). By letting x = pi = Kk\, 

y =p 2 = hk 2 , H x = M Pl = h^k^ = M Pl = cr = a p = ha k in (J72) and assigning fx pi ->■ (p 1 ) o = h(ki} = hk , 

fj, p2 — ¥ (p 2 ) Q = ^(^2)0 = —hk , o p — y (7 — hako-, we can identify the non-correlated probability distribution (|72[) with 



before , 



(fci,fc 2 ) due to That is, 



pbeforc p 



non-corr. • 



The information geometry associated with the two-particle system before collision is specified by metric (|73[) . 

2 



In a similar manner, the probability density 
distribution Comparison of (|4"2"|) and 



V> after (fci,fc 2 ,;) 



in (|42p is approximated with the Gaussian probability 



implies that when 



'QM 



< 1, 



Pi 



aftci 
QM 



(89) 



where 
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with 

r = r QM - (91) 

The expression of tqm in terms of physical quantities is given by (|4Tj) . The information geometry associated with the 
two-particle system after collision is specified by metric (|70j) . 



B. Purity as a Measure of Quantum Entanglement 

The subsystem purity of a composite system of two particles engaged in a head-on collision was calculated in [24] by 
deriving the two-particle wave function modified by s-wave scattering amplitudes. They utilized the purity function 
V as a measure of entanglement. Formally, the purity function is defined as 

V = Tr(pl), (92) 

where p 1 = Tr 2 {P12) is the reduced density matrix of particle 1 and p 12 is the two-particle density matrix associated 
with the post-collisional two-particle wave function, given by (|19[) . For pure two-particle states, the smaller the value 
of V the higher the entanglement. That is, the loss of purity provides an indicator of the degree of entanglement. 
Hence, a disentangled product state corresponds to V = 1. We remark that although V shares similar features to 
entropy, it is more readily accessible to theoretical analysis (13 , [H, [H[ than the latter. In atomic physics, V has also 
been employed to quantify two-body correlations in a multitude of dynamical processes (l2l . Il3l . [Til . |63| . 
Given the system, it was found in [24j that the purity is specifically expressed as 



V= JJJJ ^( k i' k 2>*)V'(k3,k4,t)V ^ktM^r {^^d'ktd^d'ksd^. (93) 
By employing the same dimensional analysis as developed in Section [TTl we may effectively reduce (|93[) to 

V = ( ( ( ( i/j (h, k 2 , t) tp (k 3 , k 4 , t) 41* (h, k 4 , t) xp* (fc 3 , fca, t) dkidk 2 dk 3 dh. (94) 



Now, specifying the wave-packets in (|94|) by means of (|26|) and (|3 1 jl and performing the integral, one obtains 

V = 1 - 8 (2k 2 Q + a 2 ko ) R a s + O (a 2 s ) , (95) 

where the parameter a s is the s-wave scattering length, defined from / (k) w — a s for k <C 1. Employing the scattering 
cross section E = Aim 2 , we may express the purity in an alternative manner, namely, 

V = l K — k °' +OCE). 96 

The s-wave scattering can also be understood in terms of a scattering (interaction) potential and the scattering phase 
shift. Consider a scattering potential 

where V denotes the height (for V > 0; repulsive potential) or depth (for V < 0; attractive potential) of the potential 
and L the range of the potential. Then solving the Schrodinger equation with this potential for the scattered wave, 
we are led to [64j 

fc in cot (k in L) = fc out cot (fc out L + 9) , (98) 

with 



fcin = y/W^-V) ^ Q<x<Li (99) 

fcout = x>L, (100) 

where p and £ are the reduced mass and kinetic energy of the two-particle system in the relative coordinates, 
respectively, and k m and fc out represent the conjugate-coordinate wave vectors inside and outside the potential region, 
respectively. Equation ([98| together with (f99l) and (|100ll indicates that the scattering potential shifts the phase of the 
scattered wave by 9 at points beyond the scattering region. In the next Subsection we will make use of this idea to 
determine the scattering phase shift which is linked with the micro-correlation coefficient r in our statistical model. 
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C. Information Geometric Interpretation of Quantum Entanglement 



We utilize the results of our information dynamics given by equations (f8lj) , (|82[) , (|83| , (|84|) , (|85p and (|86|) to furnish 
an information geometric interpretation of "quantum entanglement", which is characterized by the purity given by 

1. Momentum- space Gaussian Statistical Models 

To achieve the above task, one joins two different charts of Gaussian statistical manifolds, one without correlation 
(before collision) and the other with correlation (after collision) . The two models can be represented by means of (|4"5)) 
and (|90[) with associated statistical manifolds (T751) and (|T0"|) , respectively. 

The set of geodesic curves for each model is represented by equations ([ST]) , (|82|) , (|83f (for the non-correlated model) 
and by equations (|84|) . (|85l) . (|86f (for the correlated model). The two sets are joined at the junction, r = 0: r < 
(before collision) for the non-correlated model and r > (after collision) for the correlated model. 

The set of curves given by equations (|5T|) . (jg2")h may be assigned to {(pib( r )) j (pabC?")) , (fbCr)}} while the 
set given by equations (|84|) . (|85|) . (|86| may be assigned to {(pi a (r)) , (??2a(T)) , (cTa^))}- The subscripts "i" and "2" 
denote particle 1 and particle 2, respectively; subscripts "b" and " a " denote 'before' and 'after' collision, respectively. 
Then we may write the following two sets of equations: for r < (before collision), 



(pib(r)) - Mi (r; 0) = -VpI + 2^1 tanh (A G r) , (101) 
(P2b(r)) = M 2 (r; 0) = ^P? + 2<rg tanh (A oT ) , (102) 



while for r > (after collision), 



(Pla(r)) - Mi (r; r) = - ^(1 - rj (pg + 2<rg) tanh (A D r) , (104) 
(P2.(r)) = M2 (r; r) = (pg + 2^) tanh (A G r) , (105) 

Mr)) = ,( T ;0 = y|pi^-^, (106) 

where A is given by (|57|) . Here we recognize that the momenta (pib( T )) an d (Pia( T )) asymptotically converge to 
y/Po + 2erg and —y/(l — r) (p% + 2ag) toward r = —00 and r = +00, respectively (the same is true for — (p2b(T~)} and 
— (P2<l(t)}) while (ahir)} and (cr a (r)) are identical and vanishingly small toward r = ±00. Furthermore, we observe 
that there is continuity between (pi/2b( r )) an d (pi/2a( T )) an d between (cb(T)) and (cr a (r)) at the junction, r = 
(see Figure [2]). 

2. Correlation vs. Entanglement: Connection Established via Scattering and Purity 

Intuitively, if the particles are not correlated (i.e. r = 0) after collision, then no entanglement should be present. In 
this scenario, the two particle system would not experience any loss of purity so that V = 1. Indeed, this is verified 
by (|95| . From (|96|) the case V = 1 requires a s = or S = 0, that is to say, no scattering. For low energy s-wave 
scattering, f (k) ~ —a s and 6 (k) = —ka s +O (k 2 ). Thus one can readily determine the requirement necessary to 
satisfy a s — or E = 0, namely 

9 = 0. (107) 

Equation f|10T[) implies the s-wave scattering phase shift must vanish if our system is non-correlated after collision. 

A question that now arises is how to determine the scattering phase shift in view of the fact that our statistical 
model is correlated after collision. Initially, we need to examine how correlations affect the momentum geodesic curve 
(pi/2( t ))- F° r this purpose we define the momentum-difference curve (p(t)) = \ [(p2(r)) — (pi(r))]. Comparison of 
the following two equations, which follow from (|101l) . (|102|) and (|104j) . (|105J> . 

(p(r;0)) = \ [<pab(r)) - (pi b (r))] = + 2og tanh (A q t) , (108) 

(p (r; r)) = ^ [(p 2a (r)) - (p la (r))] = ^(1 - rj {pi + 2^) tanh (A r) , (109) 
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Figure 2: Plots of (pi(t)), (p2(t)) and (<x(t)) before and after collision 



indicates that at any arbitrary time r > 



(p(r;0)) > (p(r;r)), 



(110) 



while both (|108|) and (|109|) share the functional argument A r. Condition (|110p implies that the correlation causes 
the momentum to reduce for any r > (relative to the non-correlated situation). This situation is analogous to 
the change in momentum caused by a repulsive scattering potential (see ((HSJ) and (jlOOp ). It is then reasonable to 
assume there exists some connection between the scattering potential and the correlation. Provided this connection 
is established, one should be able to determine the scattering phase shift in terms of the correlation via equations 
(|9"5)) . (j9"9")l and (jlOOp . In this way, one can ultimately establish a connection between quantum entanglement and the 
statistical micro-correlation. 

Recall that before collision (at the affine time — t q ) particles 1 and 2 are separated by a linear distance R a . Each 
particle has momenta p a and —p a , respectively and the same momentum spread a . Then from (|101j) . (|102j) and (|103j) 
we have 



To give an estimate of how large A t is, we assume our momentum-space wave-packets initially have very narrow 
widths compared to their momenta such that a jp a ~ 10 -3 , for example. Then by (|87|) we find A t ~ 7.254329369. 
Using (jllip . we find p ~ 0.999999 x ^/pl + 2a 2 Q , which is equivalent to a /p ~ 10 -3 . 

For arbitrary r > after collision, the system of particles 1 and 2, which initially carried momenta p a and — p Q , 
respectively at r = — r Q before collision, now carries the relative conjugate-momentum (p (r; r)) given by (|109p due 
to the correlation. As discussed above, through (|109p . (|108[) and (| 1 10[) . it is reasonable to expect the existence of 
a connection between the correlation and the scattering potential. With non-vanishing micro-correlation the wave- 
packets experience the effect of a repulsive potential; the magnitude of the wave vectors (or momenta) decreases 
relative to the corresponding non-correlated value. One may rewrite (f9"5)) . and (|100p as 



Po = 



(Plb (-To)) 



(P2b (-t )) = VpZ + 2cr o tanh(A T ) 



(111) 




(112) 



k r cot (k r L) — k a cot (k D L + 8 a ) , 



(113) 



with 



k Q — 



(114) 
(115) 
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where O = 9 (fc Q ) w —k a s — —p a s /h denotes the s-wave scattering phase shift, and k r and k Q represent the wave 
vectors with and without the correlation, respectively. The connection between the correlation and the scattering 
potential can be established by combining (|114j) and (|115j) . 
From (jllOp one finds that the correlation renders 



Then using ()114[) . ()1 15[) and fjll 6|) . we determine the scattering potential, 

h 2 kl 



V = r£ 



2pL 



2fi' 



(116) 



(117) 



Equation (|117p clearly establishes a connection between the correlation coefficient and the scattering potential: the 
correlation coefficient is the ratio of the scattering potential to the initial relative kinetic energy of the system. From 
(|117p it is evident that our interaction potential is repulsive, i.e. V > since we consider non-negative micro- 
correlations, r € [0, 1). 

With the potential determined, one can determine the scattering phase shift by combining equations (|113[) . (|114l) . 
(fPT5j) and (fPTTl) . By solving (fTT3j) for 6 a , we find 



tan(9 D = 



k Q tan (k r L) — k r tan (k Q L) 



k r + k Q tan (fc Q L) tan (k r L) 
Substituting (|116[) into (|118p and expanding the expression in k Q L and r at the same time, one obtains 



tan 6 



~(k L) 3 + -^(k Lf 



O 



(k a L) 7 



— (k Lf + O [(k L) 7 



r 2 + (r 6 ) . 



For low energy s-wave scattering, k Q L — p L/H ^ 1, one may reduce (|119p to 



tan 6 w 9 



r(k LY 



By means of (| 11 7[) and (|120p we can express the scattering phase shift in terms of the scattering potential 

2fiVk L 3 2fiVp L 3 



3ft 2 



3ft 3 



(118) 



(119) 



(120) 



(121) 



which is in agreement with [65j . 

As the scattering potential has been determined, so too can the scattering amplitude be determined. To this end, 
we write 



f(ko) 



(122) 



k Q k Q 

for low energy s-wave scattering, k Q L = p L/h <C 1. Then the squared modulus of (|122l) . by means of (|121l) . reads 

'6 



Thus, we finally obtain the scattering cross section: 



V A \f(l>\\*~ 4 ^ 2fc o L6 16V^ 2 L 6 2 

S = 4tt|/(A; )| « = — «47ra s . 



(123) 



(124) 



Equations (|93|) and above demonstrate how the entanglement can be measured from the loss of purity by use of 
the scattering length or cross section. By combining (|95|) and (| 123[) we find the purity 



8rfc Q 2 (2fcg + al) R Q L 3 _ 16 M y (2fcg + <rg ) fl Q L 3 
3 3ft 2 



(125) 
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The correlation coefficient r can now be expressed in terms of the physical quantities such as the scattering potential, 
the scattering cross section and the purity. Solving equations (|117[) . (|124p and (|125[) for r, we obtain 



V 



2/jV 2/j,V 



£ h 2 k 2 

3(1-7?) 
K 8kl{2kl + al )R Lf 

In view of (gTJl, ([9"T]). (|124p and (|127p . one obtains the following relation: 

V _ 4ft 2 fc Q 4 (2fc 2 + a 2 ko ) R 



L 3 



3/i 



(126) 
(127) 
(128) 

(129) 



which indicates that the uniform scattering potential density is solely determined by the initial conditions of the given 
system. 

From ([I08]) . {109| and (fTT0| it is observed that for the micro-correlated Gaussian system considered here, more time 
is required to attain the same momentum value compared with the non-correlated Gaussian system. For example, 
in order to attain the same value as the initial momentum p a , the non-correlated system and the micro-correlated 
system would require time intervals r Q and r*, respectively, where 



Po 



y / p 2 + 2<T 2 t&nh(A T ), 



Po = y/(l - r) (pg + 2a 2 Q ) tanh (A t„) . 

Combining (|130l) and (I13ip . we obtain 

tanh (A q t*) = (1 — r)^ 1 ^ 2 tanh (A t ) . 
Rewriting and expanding both sides of (|132p . we have 



1 - 2e- 2A ° T * + O (e- 4A ° 
Rounding f| 133[) off and arranging terms, 



(1 



-1/2 



[1 - 2e 



-2A a 



(1 



-1/2 



.-1/2 



„2A 



-4A r 



)] 



-1/2 



The first term on the right hand side of (|134[) can be approximated to 1 since (1 - 
r <C 1. However, r in the second term should not be disregarded in the same way because (1 — r 
\\r + O (r 2 )] e 2AoTo is not negligible. Therefore, we may rewrite (|134j) as 



(130) 
(131) 

(132) 

(133) 

(134) 
O (r 2 ) and 



- 1 



o2A Q 



(1-r) 



-1/2 



(135) 



where A 

2 

exp 



r represents a new quantity that we term "prolongation", and rj A 



te) ; 



l p 2A r 
2 L 



Po 



/A, 



for ^ < 1 due to ^HTJ) . From ^35]) we find 



A oc 



lnU 



(1 



-1/2 



(136) 



At this juncture we emphasize the following points: 



• The upper bound value of r depends on the initial conditions, namely p Q and a a through the right-hand side of 
(|135p . The right-hand side of (|135l) must always be positive, so that given r -C 1, we require 



r < 



'/A 



exp 



O 



(137) 



for ^jcl. For example, with a /p ~ 10 3 , for the right-hand side to be positive we must have r < 2 x 10 6 . 
In view of (|125p . equation (|137p provides a lower bound estimate of the purity for a system with well- localized 
wave-packets, i.e. ^ < 1. 
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• With r being close to the upper bound value, A would be infinitely large due to (|136l) . On the other hand, with 
r vanishing, i.e. no correlation, A would vanish. This implies that A may serve as an indicator of quantum 
entanglement. 

• With r held fixed, A depends on the initial conditions p and er through (| 135[) . 

From the above points one may infer that the prolongation A could represent the duration of quantum entanglement 
for a given micro-correlated system and that further, the duration can be controlled by the initial conditions p a , a Q 
and the micro-correlation coefficient r. From (|135[) it is anticipated that the maximum duration would be obtained 
when r is the greatest, i.e. the micro-correlation is the strongest and the ratio a /p is the smallest (see Figure [3]). 
We emphasize that the prolongation serves to quantify the time required by a micro-correlated system - relative to a 



Figure 3: A = r* — r for different values of r: ri(the first plot) < r2(the second plot) 

corresponding non-correlated one - to attain the same momentum value (relative to the same initial reference time). 
The occurrence of a non- vanishing prolongation is in fact due to the existence of micro-correlations and therefore, due 
to the existence of scattering phase shifts. In other words, in the absence of scattering there is no time difference. 
This can be stated in yet another way as follows: "The prolongation encodes information about how long it would take 
an entangled system to overcome the momentum gap (relative to a corresponding non-entangled system) generated 
by the scattering phase shift. The entangled system only attains the full value of momentum (i.e. the momentum 
value as seen in the corresponding non-entangled system) when the scattering phase shift vanishes. For this reason, 
the prolongation represents the temporal duration over which the entanglement is active". 

From (|126p . (|127p and (|128p we observe that the micro-correlation coefficient r is directly associated with the 
quantum scattering process, and thus with the quantum entanglement. For example, the cross term (pip 2 ) in the 
definition of the micro-correlation coefficient r may represent the average interference between transmitted/reflected 
modes in the momentum degrees of freedom of particles 1 and 2. This may be viewed from a different perspective 
when considering the definition of micro-correlations (|138[) . For our statistical system in which a Pl = a P2 = cr, the 
micro-correlation coefficient reads 

r = r {pi ,p 2 ) ^ ^ ZM <*»> with a = ^((p 1/2 -( Pl/2 )Y). (138) 

Here the numerator is defined as covariance 

Cov (pi,_p 2 ) = (piPi) ~ (pi) ip2) , (139) 

and does not vanish if the statistical system is micro-correlated. In other words, if our statistical system models a 
quantum scattering process, then the relevant physical information such as scattering potential and scattering cross 
section should be encoded in Cov (p±, p 2 ). 
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With r<l, we may split our micro-correlated information geometry (|70l) into two pieces, 

, / 1 + r 2 + O (r 3 ) -r + O (r 3 ) \ 
9ab = -a -r + O (r 3 ) 1 + r 2 + (r 3 ) - 5o ab + K b , (140) 
CT \ 4/ 



where 



1 
1 
4 



9oab = ^ | I | . (Ill) 



and 



1 

Tr2 




-r + 0(r 3 ) r 2 +0(r 3 ) |. (142) 




This decomposition of the micro-correlated information geometry may provide a different (inherently IG) perspective 
on the phenomenon of quantum entanglement. In this view, the non-correlated geometry (|14ip is perturbed due to the 
presence of the quantum scattering, the information of which is encoded in (|142[) . Thus, the quantum entanglement 
manifests as this information geometric perturbation of the statistical space. 

V. CHAOTICITY, INFORMATION GEOMETRIC COMPLEXITY AND ENTROPY 

In this Section we obtain exact expressions for the information geometric analogue of standard indicators of chaos 
such as sectional curvatures, Jacobi field intensities and Lyapunov exponents. Finally, we present an analytical 
estimate of the information geometric entropy (IGE). This will lead us to uncover connections between quantum 
entanglement and the complexity of informational geodesic flows in a quantitative manner. 

A. Chaoticity 

1. Curvatures of the Statistical Manifold M^ rr , 
The Riemann curvature tensor 7Z a bcd of the statistical manifold -Mcmi 1S defined in the usual manner as [62] 

K\ cd = d c r a bd - d d r a bc + r%r f bd - r%r{ c , (143) 

where the non- vanishing connection coefficients are given in (|77p. The non- vanishing components of the Riemann 
tensor read 

1 1 r 1 

^•1212 — , A , n 7T, 7^-1313 — a 1 9 TV) 7?4323 — TT~9 TTi 72-2323 — a , 9 7T- (144) 

4er 4 [r 2 — 1) er 4 (H — 1) cr 4 (r 2 — 1) cr 4 (r 2 - 1) 

The Ricci curvature tensor lZ ab of the manifold A^ 3 °. r is defined as 

K ab = d c v ab - d b v c ac + r^ h rf d - rf c rg d . (145) 

The non-vanishing components of the Ricci tensor read 

1 r 12 

^i= 2ct 2 {r 2_ iy ^ 12 = - 2 a 2 (r 2 -l) =U2U 7122 = 2a 2 (r 2 - 1)' ^ 33 = (14(il 

Finally, we compute the Ricci scalar curvature 1Z of the manifold Ai 3 ^. 



K = K a a = K ab g ab = ~=J2 K ^ e i) ■ ( 14? ) 

1 i±J 
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That is, the scalar curvature is the sum of all sectional curvatures K (ej, ej) of planes spanned by pairs of orthonormal 
basis elements {e a = <9i?<*(p)} of the tangent space T p M^ rl . with p e M^ rr [62ll66l|. 

K M sv (u,v) = nabcdU v u v — — ; U ^rh\ v^h? with i 7^ j, (148) 

(.9arf56c - 5ac56d ) U a V»U c V d 

where (e a ,/i fc ) = (5 b a . Notice that K M 3r> completely determines the curvature tensor. The components of the 
sectional curvature are given by 

= -\ = K-^, K, 2 = - i = K-n 3 , K a = -\ = K- a . (149) 

From above, it is worthwhile to note that both statistical manifolds M.^ rr and -M^ n . COII are negatively curved, 
with the micro-correlation independent Ricci scalar curvature TZj^sd = — | = TZ^d and sectional curvature 
K M 3d = —j = K M 3d . Moreover, the constancy of the sectional curvature in all directions imply that both 
M^mT an( i -^non-corr are isotropic manifolds. Below, this will be verified by the vanishing of all components of the 
Weyl projective curvature tensor W a bcd defined on each space. 



2. Anisotropy and the Weyl Projective Tensor 

The anisotropy of the manifold underlying system dynamics plays a crucial role in the mechanism of instability. In 
particular, fluctuating sectional curvatures require also that the manifold be anisotropic. The Weyl projective tensor 
quantifies such anisotropy and is defined as 

Wabcd d = TZ a bcd r {TZbd9ac - TZ bc g ad ) , (150) 

71—1 

where n is the dimension of the manifold on which W a bcd is defined. By direct computation using (|70[) . (|144[) and 
(|146l) . we find that all components of (|150l) with n — 3 are vanishing. The fact that W a bcd — implies the manifold 
A^cS-r i s isotropic. If the manifold over which a system evolves is maximally symmetric, then 

TZ 

Ti-ab = —9ab- (151) 

n 

This is obtained from (|150j) . using the fact that W a bcd — 0. By inspection of (|146|) . (|70j) and the fact that TZ = — | 
and n = 3, it is evident that (|15ip is indeed valid for our statistical manifold M^ rl . Upon substitution of (|151[) into 
(|150p we obtain 

TZ 

Wabcd = TZabcd } 7T (9bd9ac ~ 9bc9ad) ■ (152) 

n (n — 1) 

By the fact that W a bcd — 0, this leads to 

TZ 

TZ a bcd = — -, 7T {9bd9ac ~ ffbcffad) , (153) 

n{n—l) 

which again proves to be true for our manifold A4^ TI by inspecting (|144[) . (f7TJ|) and the fact that TZ = — | and n = 3. 
Contracting the both sides of (|15ip . one finds 

5\ = n, (154) 
which is also the case for our manifold Aicon > ^ na t is, n = 3. 

3. Jacobi Fields and Lyapunov Exponents 



r \ a — 

For the sake of clarity, consider the behavior of a family of neighboring geodesies { $'L,3d (t; <f) \ on the 

statistical manifold -M|!° r , where r and if = [q 1 , <? 2 , <^ 3 ) are affine parameters. The geodesies i?^43d (t; <?) are solutions 
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of equation (T75|) . The relative geodesic spread on M.^ IT is characterized by the Jacobi-Levi-Civita (JLC) equation 

Mm, 



2 ja 



Dt 2 



£> 2 J" 



where a,b,c,d — 1,2, 3, and the second order covariant derivatives U D ^ are given by [68 
D 2 J a d 2 J a ^„ dJ b dd c „„ ,d 2 i? c „„ d?? d cW c 



pa ^ 
focJ dr 2 



Dr 2 dr 2 

and the Jacobi vector field components J a are given by 



r? 



dti f d-d c 



bc.d 



r 54* = dr{T;?) 



dr dr 



—r + nr*- — —j 



bc df dr dr 



(155) 



(156) 



(157) 



J = { J a } a=1 2 3 represents how geodesies are separating. The JLC equation of geodesic deviation is a complicated 
second-order system of linear ordinary differential equations. It describes the geodesic spread on curved manifolds 

of a pair of nearby freely falling particles traveling on trajectories $ a (r) and $' a (r) = f -d a (r) + (5i? a (r). Equation 
(|155l) forms a system of three coupled ordinary differential equations linear in the components of the deviation vector 
field (|157l) but nonlinear in derivatives of the metric tensor g a b (0). It describes the linearized geodesic flow: the 
linearization ignores the relative velocity of the geodesies. When the geodesies are neighboring but their relative 
velocity is arbitrary, the corresponding geodesic deviation equation is the so-called generalized Jacobi equation (69j . 
The nonlinearity is due to the existence of velocity-dependent terms in the system. Neighboring geodesies accelerate 
relative to each other with a rate directly measured by the curvature tensor lZ a bcd- 

By means of (|153[) the second term on the light-hand side of equation (|155j) can be rewritten as 



9& M* 

i^abcd o ** o 

or or 



n[n — 1) 



PabJ" 



where ||v|| = \J g a bV a v b with v a = ^r, and P a b = g a t, — u a Ub with u a = v a / ||v| and u a u a = 1. 
orthogonality between u a and J a , we have P a bJ b = Ja and thus ()158[) is now reduced to 



(158) 



Due to the 



i^abcd r* <J r\ 

or or 



n\w\\ 



Ja- 



n{n — 1) 

Using equations ([70)) and , the squared modulus of v a is computed as 
„ 2 <9tf a <9tf fc 



(159) 



Sab' 



9t 9t 
1 



(1 - r 2 )a 2 



dr 



d^ 2 
dr 



dr 



d^ 2 
dr 



da 



dr 2 



da 



Upon substitution of (|8"6")) into (|160|) we find 



= AAi 



(160) 



(161) 



where A Q is given by (|87j) . By combining (|155| . (|159|1 and (|161[) we finally simplify the JLC equation to the following 
form: 



2 TO 



£ 2 J' 
Dr 2 



QJ a = 0, 



(162) 



where 



(163) 
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which has been computed with n = 3, 1Z = — | and ||v|| 2 from (|16ip . 
The Jacobi vector field intensity is given by 

,. I : 
JM|J> r = ||J|| = (g a bJ a J ) 2 = (J a J a )'- 



(164) 



For applications of the asymptotic temporal behavior of Jm?d (t) as a reliable indicator of chaoticity, we refer to 

J M 3D 



our previous articles in references [28-30]. Defining the operator Cl = D 2 /Dt 2 and observing its action on 



leads to conclude 



^M3D r = (^°) Ja + ^ faa) = -2QJ a J a = -2QJ 2 M , ?iv 



(165) 



Equation f| 165|) follows from (|162p and (|164p . We may however, write ^•J^n — 2J M 3d {sijj^zv . This fact 
together with (|165p enables the further reduction of (|162p to a scalar form 



B 2 J j 



M'i 



Dt 2 



Since Q < 0, the solutions of equation (|166|) assume the form 

^M%° tr . ( r ) = ~JTq U (°) sinh (v 7 -*? 1 ") > 

dJ"^3D (r) 



(166) 



(167) 



where w (0) 



T = 



Recalling the definition of the hyperbolic sine function, sinhx = | (e 1 — e a= ), it is clear that the geodesic deviation 
on jM^St 1S described by means of an exponentially divergent Jacobi vector field intensity Jm 3D i a classical feature 
of chaos. In this Riemannian geometric approach, the quantity A»3d defined as [47l | 



\m 3D = f nm — m 

corr. r _^. 00 j- 



\j M 3D (t)| + 



dJ M 3D (r) 



|^x 3D (0)1 + 



d >7M3D r (t) 










T = 



(168) 



would play the role of the conventional Lyapunov exponents. In order to evaluate (|168p we use (|167p to find 

I^MJS, ( T )| 2 = U -q sinh 2 (%/— Qt) and |rfi7,M|° r ( r ) /^ r | 2 — ^ 2 (0) cosh 2 (y/— Qr). Thus, for the case being consid- 
ered, 



lim — In 

T — ^ oc 7" 



= lim — In 

T — ► OG 7" 



— — sinh 2 

(i-Q) ■ 



Qr ) + cosh 

2 v /r Q- 



-Qr 



Therefore, it follows that 



X M 3u T =°° 2y^Q = 2A Q > 0, 



(169) 



(170) 



which is due to (|163p . From (|170p we observe the following points: the information about chaoticity encoded in the 
positive Lyapunov exponent does not depend on the statistical correlation, i.e. = X M 3o = X M 3o = 2A a , 



and the Lyapunov exponents can be determined solely from the initial conditions (see equation (|87p ). 

B. Information Geometric Complexity and Entropy 

We recall that a suitable indicator of temporal complexity within the IGAC framework is provided by the information 
geometric entropy (IGE) S M 3d (r) [45l. |46|. 



del' 



S M 3n (t) = lim \nV M 3 



^(geodesic) ^ 



(171) 
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The information geometric complexity (IGC) is defined as the temporal average of the dynamical statistical volume 



v (geo a e S1 c, (r) def ^ / 1 , dT , vd v ^ coacslc) (t/) J (1?2) 



r—yoo \ T 







The extended volume vol 
by 

vol 



^(geodesic) o ^ ^ g effective parameter space explored by the system at time t' is given 



2^ V) = / / / P{M ™ a) (V\$\d z )dd l M 2 dd\ (173) 

/ I /^(geodesic)^ corr..H; 



where p^3D g j (z? 1 , $ 2 , $ 3 ) is the so-called Fisher density and equals the square root of the determinant g = |det {g a b)\ 
of the metric tensor g (-i? 1 , d 2 , $ 3 ) , 

P(*ffi,.*> (<?W 3 ) = yfi(?^&). (174) 
The set 2)te eodeslc ) represents a subspace of the whole (permitted) parameter space T>% otaV) in ([50|), 



^(geodesic) (r/) = |Q = (^1^3^) ; r (0) < I?" < r (r')} , (175) 

where a = 1, 2, 3, and tf a = tf Q (s) with < s < t' such that z? Q (s) satisfies ([75]). The elements of x>(s° odosic ) ( r ') are 
the macrovariables {©} whose components i9 a are bounded by specified limits of integration d a (0) and d a (r'). The 
limits of integration are obtained via integration of the set of coupled nonlinear second order ordinary differential 
equations characterizing the geodesic equations. In the case of the statistical manifold of three-dimensional Gaussian 
probability distributions parametrized in terms of = (/z^/^, <x), the integration space £)te eodeslc ) ( r ') m (| 1 73[) is 
the direct product of the parameter subspaces T ll2 and Z CT , where in the Gaussian case, Z Mi = (— 00, +00) , 
T^ 2 = (—00, +00)^ and T a = (0, +oo) CT such that 

^(geodesic) = I ^ q x ^ % X„= [(-OO, +00) ® (-OO, +Oo) ® (0, +0O)] . (176) 

In the IGAC, we are interested in a probabilistic description of the evolution of a given system in terms of its 
corresponding probability distribution on M.^ IT which is homeomorphic to 2?(s codeslc ) _ We are interested in the 
evolution of the system from Tinitiai = to rg na i = r. Within the probabilistic description, investigating the evolution 
of the system from Tj n jti a i = to Tfi na i — r is equivalent to studying the shortest path (or, in terms of the ME method 
(52l [53ll57l[58|. the maximally probable path) leading from (0) to (r). 

Formally, the IGE Sj^sd (t) is defined in terms of an averaged parametric 3-fold integral (r is the parameter) over 
the three-dimensional geodesic paths connecting (0) to (r). In the present IG approach, the IGC represents a sta- 
tistical measure of complexity of the macroscopic path = f (r) on M^mr connecting initial and final macrostates ©i 
and 0p, respectively. The path (r) is obtained via integration of the geodesic equation on M^mi generated by the 
universal ME updating method. At a discrete level, the path (r) can be described in terms of an infinite continuous 
sequence of intermediate macroscopic states, (r) — [0i, 0g_i, 0^., ©fc+i, 0f] with 0j = (tj), determined 
via the logarithmic relative entropy maximization procedure subjected to appropriately-specified normalization and 
information constraints. The nature of such constraints defines the (correlational) structure of the underlying proba- 
bility distribution on the particular curved statistical manifold M.^ ZI . In other words, the correlational structure that 
emerges in our IG statistical models originates in the information pertaining to the microscopic degrees of freedom 
of the actual physical systems. It is finally quantified in terms of the intuitive notion of volume growth via the IGC 
or alternatively in entropic terms by the IGE. The IGC is then interpreted as the temporally averaged volume of the 
statistical macrospace explored by the system, in the asymptotic limit, in its evolution from 0i to 0f. Otherwise, 
upon a suitable normalization procedure that makes the IGC an adimensional quantity, it represents the number 
of accessible macrostates (with coordinates living in the accessible parameter space 

^(geodesic) (r)) explored by the 

system in its evolution from ©i to Op. 

The temporal average in f| 1 72[) has been introduced in order to smear out the possibly very complex fine details 
of the entropic dynamical description of the system on A^St ■ Thus, we provide a coarse-grained-like inferential 
description of the system's chaotic dynamics. The long-term asymptotic temporal behavior is adopted in order to 
properly characterize dynamical indicators of chaoticity (for instance, Lyapunov exponents, entropies, etc.) eliminating 
transient effects which enters the computation of the expected value of (|173[) . In chaotic transients, one observes that 
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typical initial conditions behave in an apparently chaotic manner for a possibly long time, but then asymptotically 
approach a non-chaotic attractor in a rapid fashion. 
In the case under investigation, f| 1 72[) is given by 



a*i (">"') M 2 ( r ') cr ( 1 "') 



V 



V, 



(geodesic) 

e 



(r;r) 



= lim — 

r— >oo T 



dr' 



(177) 



M°) M 2 (0) CT (°) 



where the geodesic paths (r) = (/ij (r; r) , /i 2 (r; r) , <r (r; r)) are given in 
the information metric reads 



and 



and the determinant of 



g = det {g ab ) 



(1-r 2 ) a 6 



Substituting (|178|) into (|177p , and evaluating the integral by means of 
correlated Gaussian statistical models: 



(178) 

and we obtain the IGC for the 



V 



^(geodesic) (r . r) 



1 



lim 



1 



dr' 



Mr")i;»=o M2(r")i;»=o 



a 2 (t") 



t"=0 



' 1-r r 
1+r 



A4 3D 



3 1 sinh (X M 3vt) tanh (iA^voD-r) 
--A M 3n + 1 

4 4 T T 



(179) 



where A Q has been replaced with ^X M 3o due to (|170|) . For non-correlated Gaussian statistical models the IGC 
becomes 



V 



Ml 



^(geodesic) (r . 0) 



3 1 sinh (A_v[3dt) tanh (^A^Dr) 

--A M 3D + + 



(180) 



Inserting (|179p into (|17ip and working through some calculations, we obtain the IGE for the correlated Gaussian 
statistical models: 



1-r 
1 + r 



For non-correlated Gaussian statistical models the IGE becomes 

^M 3D ( T ;0) = AjV^dt — In (\ M 3vt) 



(181) 



(182) 



By means of (I179[) and (|180[) we compare the asymptotic (long-time limit) expressions of the IGCs in the presence 
and absence of micro-correlations, respectively, to obtain 





^(geodesic) (r . r) " 






pfeeodesie) ^. Q) 



1-r 
1 + r' 



(183) 



From (fWj) and (fT52"]) we also find 



5 » 3 s. r . ( t ; r ) - 5 A^?°„_ corr . ( r ; °) = 2 ln 



1-r 
1 + r 



(184) 



From (|183[) and (]184[) we find that both the IGC and the IGE decrease in presence of micro-correlations. In particular, 

the IGC decreases by the factor yj < 1 for r > whereas the IGE decreases by | hi fi+r) < for r > 0. 

It is evident from (| 184)) that in presence of micro-correlations the IGE is attenuated in a correlation-dependent 
manner: Sj^3d (r; r) decreases as the magnitude of the correlation increases. It is important to observe that this has 

no relation to the asymptotic (long-time limit) feature of the IGE. The correlated IGE is reduced by \ ln (■jqpf ) < 
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for r > 0, which is independent of the evolution of the system (see equation (|184[l ). When the micro-correlations 
vanish (i.e. r = 0), we obtain the expected result S M 3r> (r; r) = S M 3d (t; 0). 

..(geodesic) 



With V M 3D 



where 



(r; r) in hand, we make the following observations. From f| 183[) we find 

AC 2 

'-'total 



AC = < V 



and 



- 1, .hii | ^ 
Combining (|125p and (|185p . we obtain 



V 



(geodesic) 

e 



V, 



(geodesic) 

e 



(r;0)]y 
(r;0)]y 



- V 



v, 



(geodesic) 

e 



V 



25, 



(geodesic) 

e 



(r;r)]Y 



AC 2 

'-'total 



(185) 



(186) 



(187) 



(188) 



where the dimensionless coefficient r\ c = ^k 2 (2k 2 + cr^ ) R L 3 . From (|188[) it is evident that quantum entanglement 
and the information geometric complexity are connected. It turns out that when purity goes to unity, the difference 
between the correlated and non-correlated information geometric complexities approaches zero. 



VI. FINAL REMARKS 



In this article, micro-correlated and non-correlated Gaussian statistical models were used to model the entanglement 
of a quantum mechanical system generated by an s-wave scattering event. The IGAC was used to analyze our specific 
two- variable micro-correlated Gaussian statistical model. The manifolds M^ rl and A4^ n . COTI were used to model the 
quantum entanglement induced by head-on elastic scattering of two spinless, structureless, non-relativistic particles, 
each represented by minimum uncertainty wave-packets. The degree of entanglement was quantified by the purity 
of the system. The purity V for s-wave scattering was found in terms of the micro-correlation coefficient r, the 
interaction potential range L, the initial separation R between particles, the initial momentum p Q — hk Q and initial 
momentum spread a Q = fujko- The scattering phase shift 9 as well as the scattering cross section E were both found 
to be defined in terms of r, L and p . For r — 0, 9, £ and V (interaction potential height) are each zero while V = 1 
(indicating that the system is not entangled). The micro-correlation coefficient r, a quantity that parameterizes the 
correlated microscopic degrees of freedom of the system, can be understood as the ratio of the potential to kinetic 
energy of the system. When r ^ the wave-packets experience the effect of a repulsive potential; the magnitude of 
the wave vectors (momenta) decreases relative to their corresponding non-correlated value. The upper bound value of 
r depends onp Q and a in such a manner that r increases as p decreases. This result constitutes a significant, explicit 
connection between micro-correlations (the correlation coefficient r) and physical observables (the macrovariable p Q ). 
The role played by r in the quantities V, S, 9, and V suggests that information about quantum scattering and 
therefore about quantum entanglement is encoded in the statistical micro-correlation, specifically in the covariance 

term Cov (pi,P2) = (P1P2) — (pi) (P2) appearing in the definition (j!38[) of r. 

In summary, we proposed that the emergence of scattering-induced quantum entanglement can be understood by 
considering pre and post-collisional quantum dynamical scenarios as macroscopic manifestations emerging from appro- 
priately chosen statistical microstructures. In this view, the information geometry associated with the post-collisional 
statistical microstructure can be modelled in terms of a weak perturbation of the information geometry relative to the 
pre-collisional microstructure. In particular, quantum entanglement may be interpreted as a perturbation of statistical 
space geometry: the non-correlated geometry (|141l) is perturbed due to the presence of the quantum scattering, the 
information of which is encoded in the statistical micro-correlation terms present in (|142p . Indeed, in the case where 
r = 0, the perturbation matrix (I142p is null. Thus, the quantum entanglement manifests as a geometric perturbation 
of the statistical space in analogy to the interpretation of a static gravitational field as a perturbation of flat space. 
The perturbation of statistical geometry occurs in the 2D momentum subspace spanned by basis vectors e% = cL and 
C2 = d^ 2 . In particular, after scattering the two particles maintain a correlation among their microscopic momentum 
degrees of freedom regardless of the extent of their separation in statistical space. This fact, together with the time- 
independence of the statistical geometry [i.e. the information metric is Riemannian (rather than pseudo-Riemannian) 
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since its signature is positive definite (rather than positive semi-definite)] leads to a notion of statistical non-locality. 
The perturbation of statistical geometry is associated with the scattering phase shift in the statistical momentum 
space. 

The prolongation, denoted A, was defined as the time required for the observed momentum difference between 
a correlated and corresponding non-correlated system to vanish. The prolongation encodes information about how 
long it would take an entangled system to overcome the momentum gap generated by the scattering phase shift. The 
entangled system only attains the full value of momentum (i.e. the momentum value as seen in the corresponding 
non-correlated system) when the scattering phase shift vanishes. For this reason, the prolongation represents the 
temporal duration over which the entanglement is active. It was found that for r values close to its upper bound, 
the prolongation A becomes infinitely large. On the other hand, with r vanishing (i.e., no micro-correlation) A is 
identically zero. With r fixed however, the prolongation A depends on p and a Q . Thus, the prolongation A may 
be taken to represent the duration of quantum entanglement for a given correlated system where the entanglement 
duration can be controlled by the initial conditions p a and a Q as well as r. Maximal prolongation occurs when r is 
greatest and the ratio a /p is smallest. For small initial r and p , A would be correspondingly small, suggesting that 
for such scenarios quantum entanglement is transient. 

It was determined that both statistical manifolds M^ rt and ■M^ a _ COTI are negatively curved, with a micro- 
correlation independent Ricci scalar curvature IZ^u = — | = TZj^sd . Moreover, the sectional curvature 
throughout both manifolds was determined to be constant, K M 3o = — | = K M 3r> . The constancy of the 
sectional curvature in all directions imply that both M^ TI and ■M^ n . COI1 are isotropic manifolds. This was verified 
by the vanishing of all components of the Weyl projective curvature tensor W a bcd defined on each space. The com- 
plexity of geodesic paths on A4^ rr and Al^o n _ corr was characterized through the asymptotic computation of the IGE 
and the Lyapunov exponents on each manifold. The Lyapunov exponents in both cases were found to be constant and 
positive definite, i.e. Xj^ao = Xj^ao — 2A > 0. The IGE Sj^sd (t; r) in presence of micro-correlations assumes 
smaller values relative to the non-correlated case Sj^3n (r; 0) while the growth characteristics of both correlated 
and non-correlated IGEs were found to be the same. Specifically, the larger the micro-correlation (i.e. the closer r is to 
1) the lower the values of the IGE S M 3o (r; r). Thus, the stronger the micro-correlation, the larger the gap between 
Sj^au (t; r) and Sj^au ( T ;0). This implies that £^30 (t; r) < S_^3d ( t ;0). When micro-correlations vanish 
(i.e. when r = 0), we obtain the expected result, (t;0) = r ( r ;0). In the model investigated in this 

work, the appearance of micro-correlation terms in the elements in the Fisher-Rao information metric leads to the 
compression of V^3D 2?(g codeslc ) r ^ the fraction »J -j-^ and thus, to a reduction of the complexity of the path 

leading from initial macrostate Oj to final macrostate ©p. 

Information Geometry and Maximum Relative Entropy methods hold great promise for solving computational 
problems in classical and quantum physics. Our theoretical formalism allows for the analysis of physical problems 
by means of statistical inference and information geometric techniques, that is, Riemannian (differential) geometric 
techniques applied to probability theory. The macroscopic behavior of an arbitrary complex system is a consequence 
of the underlying statistical structure of the microscopic degrees of freedom of the system. We are confident that 
the present work represents significant progress toward the goal of understanding the relationship between statistical 
micro-correlations and quantum entanglement on the one hand and the effect of micro-correlations on the dynamical 
complexity of informational geodesic flows on the other. It is our hope to build upon the techniques employed in this 
work to ultimately establish a sound information geometric interpretation of quantum entanglement. 
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Appendix A: Integration of the Geodesic Equations 

The coupled ODEs ([75]) . (|?TB) and ([5D1 can be solved via the following strategy. First, (|75j) and ([79)) can be rewritten 

as 

Ml/2 _ 2o^ 



(Al) 
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where " denotes a differentiation with respect to r. (|A1[) can be recasted as 

2a' 



x 
x 



Moreover, since 



we find 



-, where x = fi[/ 2 an d %' = Mi/2- 



x' 2a' d , J, 

— = => — In x = 2— In a, 

x a dr dr 



(It 



In 



dr = 2 



dr 



In a ) dr ln\x\ + k = 21n<7, where k — const. 



Exponentiating both sides of the above equation leads to 

exp (In |x| + k) = exp(21nu) e k \x\ = a 2 . 

Thus, 

|Ci/2 1 cr , 



Ml/2 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



where C1/2 are the integration constants corresponding to Mi/2- In order for our Gaussian statistical model to have 
smooth and natural evolution, a(r) must be positive definite and well-behaved (continuous and differentiable) over 
the entire domain of r; r 6 (—00, +00). Then from (|A6|) M1/2 mus t be either positive definite or negative definite over 
the entire domain of r and the sign of Ci/ 2 must be associated with the sign of /i'1/2 so that a is positive definite and 
free from nodes. We can rewrite (IA6I) as 



Substituting (|A6|) into (|80| . we obtain 



M1/2 = C 1/2 a 2 \ with — ^ > 0. 

W/2 



(2rCiC 2 - C 2 - C 2 ) a 3 



Dividing both sides of (| A8|) by cr yields 



4 (r 2 - 1) 



1 



{2rdC 2 - C\ - Cl) a 2 = 0. 



0-2 4(r 2 -l) 

The first term of (|A9I) can be rewritten as a complete differential by means of the following identity: 



dr \ a 



(A7) 



(A8) 



(A9) 



(A10) 



The second term of (|A9[) is also a complete differential form due to (|A6[) . Then, for /i x and /i 2 , respectively, we may 
rewrite (IA9I) as 



dr \ a 



Ci 



4( r 2 _ 
C 2 



4( r 2 _ X ) 

Integrating both sides with respect to r, these become 



Ci V Cr 

C 2 V ^2 



- 1 
- 1 



Ml = 



/'2 



a 4(r 2 - 1) 
a 4(r 2 - 1) 



Co 



I 2t — 



c 2 



C 2 



Mi + Di 

A*2 + Aj 



0. 
0. 

0, 

0, 



(All) 
(A12) 

(A13) 
(A14) 
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where D\ and D 2 are integration constants. 

Substituting (|A13j) and (|A14|) into (|A"Tj) leads to 



Mi + 

/4' + 



2(r 2 - 1) 

C 2 
2(r 2 - 1) 



C2 
Ci 

c 2 



2r 
2r 



C2 
Ci 

c 2 



^!^ + 2£>i/ii = 0, 
M2M2 + 2D 2 /4 = °- 



Then integration of both sides of (|A15|) and (|A16[) with respect to t yields 



/'2 



1 4(r 2 - 1) 
C 2 

4 ( r 2 _ ]_) 



"C 2 

Cl 


^2r- 


c 2 \ 
cj 


- 1 


Mi 


c 2 


^2r- 


c 2 y 


- 1 


M2 



/if + 2D 1 fi 1 + E 1 = 0, 
nl + 2D 2 n 2 + E 2 = 0. 



Equations (|A17p and (|A18|) can now be represented by the general form: 

[i + a/j 2 + b/i + c = 0, 



(A15) 
(A16) 

(A17) 
(A18) 

(A19) 



which is known as the "Riccati equation" [4 
may be modified to a more tractable form: 



where 



Due to the fact that a, b and c are all constants in our problem, (|A19|) 



v' + Au 2 + B = 0, 



;/ = [i - 
A = a, 



2a 



The solution of (|A20[) is given by the form: 



B = -—+c. 
4a 



1 au' + f3v' 



A au + /3v ' 

with a and /? being arbitrary constants, not both zero, while u and v are linearly independent solutions of 



d ( 1 dz\ 
aW \~A~aW) 



d 2 z 

}+Bz = => —r + ABz = 0. 
dr z 



One finds easily 



where 



u 

V 



7 = ±\/-AB = ± 



\Jb 2 - 4ac 



(A20) 

(A21) 
(A22) 

(A23) 
(A24) 
(A25) 



(A26) 
(A27) 



(A28) 



Then by means of equations (|A2i~j) . ([A22]) . ([A23|) . (|A24| . (fA26|) . (|A27|) and (|A28ll . we find 



M(r) 



7 ae 



7r 



--JT 



a ae^ T + f3e-~< T 2a 
± 



Vfe 2 -4ac aex P 


1 Vfc 2 -4ac_ 
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— /3exp 
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=F V 2 T J 


^ a a exp 


1 Vfc 2 -4ac_ 
- 1 - 2 ' 


+ /3exp 


-,- Vb 2 -4ac_ 
2 ' 



b 

2a' 



(A29) 
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Finally, we may identify equations (|A17|) and (|A18I) with (|A19|) to find the solutions fi ± and (i 2 via (|A29I) : 

7i/ 2 a 1/2 e^ T - P 1/2 e-^ T b 1/2 



M1/2 ( T ) 



where for ^ 



with 



1/2 ai/ 2 e 7l / 2T +P 1/2 e 7l /^ T 2a 1/2 

_ , \Jb\ - 4a!d 
7i = ± 7, , 



a x = 



4 (r 2 - 1) 
bi = 2D U 



C 2 



Co 



2r - - 1 



Ci 



and for // 2 



with 



a 2 



b 2 = 



7 2 = ± 



C 2 
4( r 2 _ 1) 



_ , V & 2 ~ 4a 2 c 2 



Ci 

c 2 



2r 



Ci 

c 2 



C2 = #2 



(A30) 



(A31) 



(A32) 

(A33) 
(A34) 



(A35) 



(A36) 

(A37) 
(A38) 



In order for our system to have non-oscillatory and non-constant evolution, 7]y 2 must be real, thus the quantities 
b\, 2 — 4ai/2Ci/2 must be positive definite. Later, we will find the conditions for this (see (|A53|) ) . 
We may rewrite (|A30[) as 



Ml/2 0) = 



7i/ 2 <5 1/2 e 7l /2 r - e 7 */2 



bi/2 



a 1/2 8 1/2 e^/* T + e-^/* T 2a 1/2 



where 



/2 



H/2 



Pi 



/2 



(A39) 



(A40) 



By means of (|A7p and (|A39|) we find 



a(r) 



'Mi/2 ( T ) 



Ci 



/2 



<5l/2 



0.1/2C1/2 



7l/2 



(A41) 



where ai/2, 71/2, ^1/2 ar e given by (| A32|) . (|A36j) . (|A31ll . (|A35|) and (|A40|) . However, our a obtained either via n 1 or 
via n 2 must be identical. This yields the following equality: 



A8 2l 2 2 



Ci C 2 aiCi(5ie~<i T + e~~<i T ) 2 a 2 C 2 (5 2 e^ T + e~^) 2 ' 



In order for (jA42[) to be generally true, the following conditions must be satisfied: 

iTxl = |7 2 |, 
Si = s 2 = 1, 

aiCi — a 2 C 2 . 



(A42) 



(A43) 
(A44) 
(A45) 
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From (|A32[) and (|A36[) one finds that (|A45I) holds true by itself. In order for (IA43[) to hold true, we require 



D 2 CtEh 



i 4(r 2 -l) 



9i 
c\ 



2r- 



9i 

Ci 



- 1 



= Di- 



C2E2 
4(r 2 - 1) 



C 2 \ C 2 



- 1 



Substituting the conditions (|A43|) and (|A44|) into (|A39|) . we obtain 

V1/2 = Z — tanh (77-) - 



Ol/2 



0.1/2 



where 7 = I7J = |7 2 |. Adding [i 1 and fi 2 , we find 



Mi + M2 



1 1 



H 7tanh(7r) h 



Di , D2 



ai CL2 



(A46) 



(A47) 



(A48) 



We make use of this Gaussian system to model a head-on collision between two Gaussian packets in momentum space, 
where each particle carries the average momentum, (pi) = [i Y an d (fa) = M2j respectively. Thus, the total momentum 
of the two-particle system represented by (|A48[) must be conserved. This requires 



1 + 1=0. 

01 a 2 



(A49) 

For convenience we require both /Xj (r) and [i 2 ( r ) cross at r = 0. From (|A47[) we find that this condition implies 

Dx = D 2 = 0. (A50) 

From (|A49p one finds 

C\ = -C 2 = C ^ 0. (A51) 
Then due to (|A45|) . (|A50|) and (|A51|I . (jA46|) is reduced to 

B a = -_E 2 e£^0. (A52) 
Substituting ([A50]) . (|A5l|l and (|A52|) into (j A46|) . we obtain the above mentioned reality condition for 7^2 ; namely 

CE < 0. 

From (lAlUl) together with (lA32l) . (|A!3l) . (|A34l . (lA43l . <|A50|) . (IA3T1) . JM2J, we find 



7 = |7xl = I72I 



CE 



, 2(r-l) 

Then substituting (|A50j) . (|A51|) and (|A54|) into (|A47p . and using (|A41[> . we finally obtain 



(A53) 
(A54) 



Mi (r;r) 

M2 far) 
cr(r;r) 



2£(Y - 1) 



c 



tanh 



CE 



2(r-l) 



2£ (r - 1) 



C 



tanh 



2(r-l) 



C 



(A55) 
(A56) 
(A57) 



where we have set C < and E > 0. In our probabilistic macroscopic approach to dynamics, these geodesic trajectories 
represent the maximum probability paths on -M^. r . 

For the non-correlated Gaussian system, we set r = in (| A55|) . (|A56[) and (IA57|) to obtain 



Mi (r;0) 
M2 (r;0) 



2£ , 
- — p tanh 
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2E 
- — tanfi 
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(A58) 
(A59) 
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Distinguishing the constants C and E for the correlated Gaussian system from those for the non-correlated Gaussian 
system, we rewrite (|A55[) . (|A56[) and (|A57[) as 



*™ = f^*"* Uw^Tf)- (A62) 



where the subscript " r " in C r and E r implies that the constants are dependent upon the correlation coefficient r of 
the given statistical manifold. 



Appendix B: Refining the Geodesic Trajectories 



In this Appendix we join two different charts of Gaussian statistical manifolds, one without correlation (before 
collision) and the other with correlation (after collision) . The set of geodesic curves for each model is represented by 
equations (|A58[) , (|A59[) , (|A60[) (for the non-correlated model) and by equations (|A61|) , (|A62[) , (IA63[) (for the correlated 
model). The two sets are joined at the junction, r = 0: r < (before collision) for the non-correlated model and 
t > (after collision) for the correlated model. 

The constants, C and E in (|A58[) , (|A59|) and (|A60p can be determined via the conditions at the initial affine time, 
— t - We assign the initial momenta and the dispersion of the wave-packets as 



Mi (-r o ;0) 
fi 2 (-r o ;0) 
<t(-t o ;0) 

Combining (jBlj) and (|B3j) . one obtains 
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-Po, 



E _ 1 2 
C ~ 2 Po 



Also, taking the ratio between a Q and p via ()B1[) and (|B3|) yields 

CTo f 
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°° \/2sinh( 



2 T< 



Upon considering large r in (|B5|1 . we find 
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1 ( Oo_ 

Po 



8 \Po 
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(Bl) 
(B2) 
(B3) 



(B4) 



(B5) 



(B6) 



Equation ()B6|) implies that r Q should be chosen sufficiently large so that the ratio a /p will be very small, while 
— remains finite. The constants C and E can be individually determined by simultaneously solving (|B4[) and 
(IB61). 
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In a similar manner, the constants C r and E r in (|A61|) . (| A62|) and (|A63|) can be determined via the conditions at 
the reversal time r Q . We assign the momenta and dispersion of the wave-packets according to 



M T o;0 = 
M 2 ( T o;r) = 
f(r ;r) = 

Combination of (1571) with (lB9l leads to 
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r cosh 



E T 



(\/#=Ty r °) 
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J 2 



C r 2(1- r) 

From (|B7[) and (|B9[) it is found that the ratio between cr^ and p' Q reads 

cr' 1 1 



P'o VW^r) sinh (^5r ) 
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(BIO) 
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From (|B11[) it is found that for large t q , 
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• sinh 1 
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(B12) 



Here again it is implied that r Q should be taken sufficiently large so that the ratio, cr' /p' can be very small while 
\J 2(r-i) remams finite. Furthermore, the constants C r and E r can be individually determined by simultaneously 
solving (fBlO]) and JEE}. 

The two sets of geodesic curves (with and without correlations) are joined at the junction r = 0. The two sets of 
geodesic curves must be continuous at the junction r = so as to ensure the collision does not assume any unphysical 
irregularity in the momentum dispersion. From (|A58D . (|A59jl . (|A60ll and (|A61jl . (|A62j) . (|A63|) it is found that this 
continuity condition is satisfied by 



E _ E r 

C Cr 

Using condition (|B13[) together with (|B1[) and (|B7|) . one may compare p with p' Q as follows, 

tanh 



(B13) 



Po n 

— = V 1 — r 

Po tanh 



CrE r 

2(r-l) T ° 



(B14) 



where e = exp (— v 7 — 2C£ , r ) and e' = exp 
also compare <r with a^, 



2C r E r , 



j^To J . In a similar manner, by way of (|B3[) and (|B9[) one may 



cosh 



cosh 



{l + (6- e ')+O[( £ -0 2 ]} 
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From (|B14[) it is observed that for sufficiently large r the ratio p' /p is not significantly influenced by how the 
functional arguments \l ~^-T n and 



» ' " y 2(r-i) T ° com P are with each other, since the quantities on the right-hand side 
of e and e' are very small (as is the difference e — e'). From (|B15|) however, the ratio cr' /a appears to be influenced 
by how those functional arguments compare with each other since the leading approximation reads 



(B16) 



— = exp 

e 




From (|B16|) . one observes that the difference y — -^p — y 2(r~i) mus t vanish in order for \J e'/e to remain finite given 
that r is sufficiently large; otherwise a non-vanishing difference could result in a sufficiently large exponent when 
multiplied by a large value of r Q - this would cause y/e' /e to grow or decay exponentially. A vanishing value of 
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implies 



From (jBT4|) . (jETB] and (jETTf it follows that 



Pa 



rPo, 



(B17) 



(B18) 
(B19) 



Equations (|B18|) and (|Bf9|) also satis fie s the condit ion (IB13)) th ro ugh (TBI)) an d ()B10p . 

Sub stitut i ng (pp l jB6 | and (fBT0|) . (fBf~2f into (fA58|) . (|A59|) . ((Alp) and (jA6lj) . (fA62]l . (fA63|) . respectively, and 
using (|B17|) . (|B18|) . (|B19|) . we may rewrite the geodesic trajectories as follows: for the non-correlated Gaussian system, 



/ii(r;0) = -VPo + 2a2tanh(A r), 
/i 2 (r;0) = VPo + 2CT o ta nh (4>r) , 

t(t;0) ' ' 



°cosh(A>T) ! 



while for the correlated Gaussian system, 



where 



Ml (r;r) = -y/(l-r) (rf + 2o*) tanh (A*t) : 

t 2 ( r ; r ) 

a(r;r) 



M2 (r;r) = ^(1 - r) (p§ + 2a 2 ) tanh (A oT ) 
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which is defined from (IB6I). 
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